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In this paper, we study the recent excess of low energy events observed by the CoGeNT collabo- 
ration, and discuss the possibility that these events originate from the elastic scattering of a light 
(ttidm ~5-10 GeV) Weakly Interacting Massive Particle (WIMP). We find that such a dark matter 
candidate may also be capable of generating the annual modulation reported by DAMA, without 
conflicting with the null results from other experiments, such as XENON 10. The regions implied 
by CoGeNT and DAMA are also near those required to produce the two observed CDMS events. 
A dark matter interpretation of the CoGeNT and DAMA observations favors a region of parameter 
space that is especially attractive within the context of Asymmetric Dark Matter models. In such 
models, the cosmological dark matter density arises from the baryon asymmetry of the universe, 
^ ' naturally leading to the expectation that ttidm ~1-10 GeV. We also discuss neutralino dark matter 

from extended supersymmetric frameworks, such as the NMSSM. Lastly, we explore the implica- 
, tions of such a dark matter candidate for indirect searches, and find very encouraging prospects for 

• experiments attempting to detect neutrino or gamma ray annihilation products. 
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^ ■ I. INTRODUCTION 

CIh' The direct and indirect detection of dark matter has recently become a topic of intense discussion and activity, with 
1^ . a number of reported signals having been interpreted as possible detections of particle dark matter. Among direct 
detection experiments, the DAMA collaboration, which has collected over 1.17 ton-years of data, reports an annual 
modulation in their event rate with 8.9a significance, which they interpret as evidence of dark matter The CDMS 
' experiment has also observed a small excess (2 events, corresponding to less than 2a significance) in their low energy 
I nuclear recoil window [2| . Furthermore, the CoGeNT collaboration has very recently announced the observation of an 
■ excess of events at low energies 0]. If the CoGeNT excess is interpreted as a detection of dark matter, this points to 
dark matter with a mass in the range of approximately jtidm ~5-10 GeV, and an elastic scattering cross section with 
nucleons of approximately cr ~ 7 x 10~^^ cm^. Remarkably, a dark matter candidate with this approximate mass and 
cross section would also be capable of producing both the two events observed by CDMS, and the modulation signal 
^ [ reported by DAMA. 

Indirect detection efforts have also recently produced a number of tantalizing, if unconfirmed, signals of dark matter. 
The PAMELA experiment has reported an excess of cosmic ray positrons between 10 and 100 GeV relative to the 
predictions of standard cosmic ray models Q . The Fermi Gamma Ray Space Telescope (FGST) collaboration has also 
reported a significantly harder spectrum of electrons (plus positrons) than had been anticipated Q. In addition, there 
are apparent excesses in radio and gamma rays originating from the inner kiloparsecs of the Milky Way: the WMAP 
;h ' haze [0], and the recently reported Fermi (gamma ray) haze 0. Lastly, it has recently been shown that the spectrum 
and angular distribution of gamma rays from the Galactic Center, as observed by the FGST, is consistent with an 
annihilating dark matter interpretation . Although each of these observations is consistent with annihilating dark 
matter, it is not yet possible to rule out less exotic astrophysical origins, such as pulsars 

In this paper, we focus on recent results from direct detection experiments, and examine the implications of the 
detections made by CoGeNT, DAMA, and CDMS. ^ We begin by discussing the compatibility of these observations 
with each other, and with the null results of other experiments; in particular the recent low threshold analysis of 
the XENONIO collaboration [ll|. We then turn to a discussion of dark matter models that are naturally able 
to accommodate the observations of CoGeNT, DAMA, and CDMS. Particularly attractive within this context are 
Asymmetric Dark Matter (ADM) models [H, [l^, in which the dark matter abundance is related to the baryon 
asymmetry of the universe. In such models, the dark matter's mass is predicted to be a factor of several times 
heavier than the proton; precisely within the region indicted by CoGeNT. We also explore supersymmetric models in 
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The dark matter candidates giving rise to direct detection signals considered here could not generate the positron excesses in Fermi and 
PAMELA, though multiple components of dark matter could simultaneously give rise to multiple signals of different type in direct and 
indirect detection experiments [lOll . 
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which these signals can be produced. In particular, we discuss neutralinos within the context of the next-to-minimal 
supersymmetric standard model (NMSSM). Lastly, we turn our attention to the prospects for indirect dark matter 
searches. We find that, while the region implied by CoGeNT is not ruled out by any indirect detection experiments, the 
outlook for dark matter's indirect detection is very encouraging in this scenario. In particular, neutrino and gamma 
ray searches for dark matter annihilation products are currently within a factor of ~2 of the sensitivity required to 
test many of the models capable of generating the signals reported by CoGeNT, DAMA and CDMS. 



II. COGENT, DAMA, CDMS AND XENON 



In this section, we explore whether a dark matter interpretation of the CoGeNT excess is consistent with the signals 
reported by DAMA and CDMS, and whether such an interpretation is consistent with the null results of XENONIO [T^ 
and the CDMS silicon analysis [l5|, [l^. In particular, we follow Refs. [l7l - [20j in studying direct detection signals in 
the low mass WIMP region, including the possible effects of channeling in the DAMA experiment [2l| . 



A. Detection Rates 



The differential rate of dark matter elastic scattering events per unit detector mass in nuclear recoil energy is given 

by 

dR ,j PPM f ,3 da 

——=,Nt / dh,vf(v,Ve)——, (1) 

where Nt is the number of target nuclei per unit mass, muM is the dark matter particle mass, poM = 0.3GcV/cm'^ 
is the local dark matter density, v is the dark matter velocity in the frame of the Earth, i)^ is the velocity of the Earth 
with respect to the galactic halo, and f{v, v'e) is the distribution function of dark matter particle velocities, which we 
take to be a standard Maxwell-Boltzmann distribution: 

/(.Vr.) = ^e-(--)V< (2) 

The Earth's speed relative to the galactic halo is = Vq -\- WorbC0S7C0s[w(t — tg)] with vq ~ vq + 12km/s, Worb = 
30km/s. C0S7 = 0.51, to = June 2nd, and lu — 27r/ycar. We set vq = 270km/s throughout most of our analysis. The 
upper limit of the velocity integration of Eq. ^ is the galactic escape velocity, 490km/s < v^sc ^ 730km/s at 90% 
C.L. [i^l, and the minimum dark matter velocity, Vmm, is 



ERtriN 

where Eji is the recoil energy, niN is the mass of the target nucleus, and is the dark matter-nucleus reduced mass. 
For a spin-independent cross section between dark matter particles and nuclei, we have (23l | 

da niN aN [fpZ + J^jA - Z)f ^ 

where is the reduced mass of the dark matter particle and nucleon (proton or neutron), ctat is the scattering cross 
section of the dark matter particle with neutrons, Z and A are the proton and atomic numbers of the nucleus, and fn^p 
are the coupling strengths of the dark matter particle to neutrons and protons respectively. The couplings fn,p are 
calculated from a coherent sum over the couplings to the quark constituents of the nucleon. We use the Woods-Saxon 
form factor, F{Eh) [Mill. 

CoGeNT is an ultra-low noise detector operating in the Soudan mine, consisting of 0.33 kg (fiducial mass) of 
germanium. Over a period of 56 days, the CoGeNT detector recorded approximately '--^lOO events (with characteristics 
corresponding to interactions in the bulk of the crystal) above their expected background with ionization energy 
between 3.2 keVee^ and their threshold of 0.4 keVee Q. The CoGeNT collaboration has pointed out that their 



^ Ionization energies are given in units of keVee, or keV electron equivalent. For a 1 keV nuclear recoil, for example, the equivalent electron 
energy (in keVee) is X IkeV, where is the quenching factor for the nuclear material composing the scintillator. For the germanium 
target used in CoGeNT, the quenching factor is approximately 0.2 to 0.3 for events in the energy range of the observed excess. For the 
materials composing the DAMA detectors, gjva ~ 0-3 and qj 0.09. 
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FIG. 1: The regions in the elastic scattering cross section (per nucleon), mass plane in which dark matter provides a good fit to 
the CoGeNT excess, compared to the region that can generate the annual modulation reported by DAMA at 90% confidence 
(darker grey regions). In this figure, we have adopted vo = 270 km/s and use two values of the galactic escape velocity: 
Vcsc = 490 km/s (left) and Vcsc = 650 km/s (right). In calculating the DAMA region, we have neglected the lowest energy 
bin (the effect of this is shown in later figures) and treated channeling as described in Ref. [5^ . If a smaller fraction of events 
are channeled in DAMA than is estimated in Ref. [2^, the DAMA region will move upward, toward the yellow regions (near 
ajv ^ (,jjj2^ which include no effects of channeling), improving its agreement with CoGeNT. Also shown is the 90% C.L. 

region in which the 2 events observed by CDMS can be produced. If the escape velocity of the galaxy is taken to be relatively 
large, this region can also approach those implied by CoGeNT and DAMA. Constraints from the null results of XENONIO and 
the CDMS silicon analysis are also shown. For the XENONIO constraint, we have used the lower estimate of the scintillation 
efficiency (at la) as described in Ref. [27| . 



events can be fit very well by a ~10 GeV dark matter particle with an elastic scattering cross section of ~ 7 x 10 
cm^ . In Fig. [U we confirm this conclusion, where we show the parameter space region in which elastically scattering 
dark matter can accommodate the CoGcNT excess at 90% confidence. In this figure, we have used vq = 270 km/s 
and Uosc — 490 km/s. Here, and throughout this paper, 90% (99%) confidence regions are defined as contours of 
=L^min "I" 4.61 (9.21), while constraints from null experiments are defined as 90% limits based on the maximum 
gap [281 method. To carry out this fit, we have assumed that the background is well described by an exponential plus 
constant, and we have required bin-by-bin that the background not exceed the amplitude of the dark matter signal. 
Without a constraint on dark matter signal to background, the entire spectrum is well fit by a pure exponential 
background. Tighter constraints on the amplitude of the background will correspond to the dark matter signal region 
shifting to larger cross sections. We fit the data in 0.05 keV-clcctron-cquivalcnt (keVec) bins from threshold at 0.4 
kcVcc to 1.8 kcVcc where the dark matter signal is negligible. Peaks in the data (consistent with a background from 
radioactive tin) at 1.1 and 1.29 keVee are fit by Gaussians of relative height 0.4 and with width consistent with the 
experimental resolution at those energies (0.0774 and 0.078 keVee respectively). We can see that for appropriate 
choices of the halo model and the fraction of channeled events in DAMA, the CoGeNT region can be consistent at 
90% C.L. with the DAMA signal and the null results XENON and CDMS-Si. Some consistency between the preferred 
region for CDMS with DAMA and CoGeNT can also be found. We now turn to discussing in detail how this occurs. 
The DAMA experiment [l| observes an annual modulation in their count rate, which can be parameterized as 

Ri = R° + Slcos[uj{t - to)]. (5) 

The subscript i in this expression denotes different energy bins. The constant term i?? is composed of both a signal 
component coming from dark matter initiated processes, and a background component arising from other sources of 
nuclear recoil: R° = b° + S°. The expressions for 5"° and S} are obtained by integrating Eq. (|4]) over a given energy 
bin. 

Channeling is a potentially important but difficult-to-predict theoretical effect which can significantly change the 
interpretation of DAMA's signal, especially when comparing this signal to the results of other direct detection cxperi- 
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ments. In a typical nuclear recoil event, only a fraction of the total energy is detected (as a combination of scintillation 
light, heat, and ionization, depending on the detector). The ratio of the observed energy to the total recoil energy is 
known as the quenching factor. For crystal scintillators, however, such as those used by the DAMA collaboration, a 
portion of the events will be "channeled" , causing most of the recoil energy in those events to be observed (effectively 
changing the quenching factor to q ~ 1). This occurs when the incident particle interacts only electromagnetically 
with the scintillator material, which can occur for certain energies and incidence angles. The importance of this effect 
for the DAMA experiment was first discussed in Ref . [2l| , and an analysis of its effect was performed by the DAMA 
collaboration in Ref. [1^ ; we refer the reader to these references for a more detailed discussion. 

Questions of how to quantitatively account for the effects of channeling in DAMA have made efforts to interpret 
their results, and to compare their results to those of other experiments, somewhat difficult. In particular, the 
compatibility of the DAMA modulation with the CoGeNT and CDMS excesses (and the compatibility of DAMA with 
the null results of XENON and other searches) depends strongly on the fraction of elastic scattering events which 
are channeled rather than quenched. An estimate of the fraction of channeled events based on simulations is given in 
Fig. 4 of Ref. [2^ . We use the following simple parameterization in our analysis [1^ : 

~ l + 1.14Sj^(keV)' ~ 1 + 0.75£;i?(keV) ' ^""^ 

It should be noted, however, that channeling at DAMA has not been experimentally verified, and studies of Nal at 
energies above the region of interest for DAMA have failed to find such an effect [30[. We also consider the DAMA 
preferred region if channeling is not present. In this case, scattering off of iodine cannot reproduce the DAMA signal 
at the low masses of interest. Scattering off of sodium is relevant, however, and is shown as a yellow region in each 
frame of Fig. [T] The cross sections required to explain DAMA without any channeling are somewhat too high to 
explain CoGeNT, but the range of masses needed to explain the two experiments is surprisingly similar. 

In Fig. [1] we show along side the CoGeNT region the parameter space that can accommodate DAMA's annual 
modulation signal. In addition to the DAMA and CoGeNT regions, we also show the region of the parameter space 
in which the 2 events observed by CDMS can be generated (consistent with the small excess recently reported (H). 
If the actual fraction of channeled events is smaller than its estimated value in Eq. [6l then the DAMA region simply 
shifts upward in cross section to accommodate the change. Thus, in Fig. 1(a), one can see by eye how modifying 
the channeled fraction to 60% of its estimated value easily moves the DAMA region into better agreement with 
CoGeNT. In Fig. 1(b) we also see how vq shifts the DAMA region to lower masses. The two constraint plots in 
Fig. 1 are considerably more favorable for a light elastically scattering WIMP interpretation of DAMA than previous 
results (3l| . The reason for this being that the allowed region of parameter space depends quite sensitively on 
assumptions about the halo model and the XENONIO scintillation efficiency. We next turn our attention to these 
details and how they impact the constraints from CDMS and XENONIO. 



B. Consistency with Null Results 



In the light WIMP window, the null results of XENONIO are among the most constraining, due to its relatively low 
threshold of 4.5 keV. Recently that threshold was lowered further to 2 keV [ll|, although no constraint was derived on 
low mass, elastically scattering dark matter in that analysis. To derive the XENONIO constraint, software efficiencies 
given in Ref. [ll| are coupled with an 5*2 twelve electron efficiency and an 5*1 detector acceptance efficiency taken 
from Ref. [s^l- At face value, the XENONIO constraint appears to completely rule out the regions of the low mass 
window favored by DAMA and CoGeNT. 

There are reasons to be skeptical of such a conclusion, however. In particular, a new measurement was recently made 
of the scintillation efficiency (the fraction of nuclear recoil energy that goes into scintillation light) of liquid xenon, 
LeS [131 • According to the new measurement, the scintillation efficiency is significantly lower than previously reported 
by the XENONIO collaboration, who take LefF = 0.19 in their analysis. A lower scintillation efficiency translates to 
a reduced sensitivity of xenon-based detectors to light dark matter. Although the central values of the new L^s 
measurement still lead to considerable tension with the region implied by CoGeNT, if we adopt values of L^ff which 
are la below its central values, we find that consistent regions can be found. This constraint is included in Fig. [2l 
along with the constraint from the CDMS silicon analysis |l^ [l^ . We show in Fig. [5J for comparison, the constraints 
from XENONIO using the old L^s = 0.19, along side those using the new measurement with and without taking into 
account the la uncertainty in each bin. Note that the XENONIO excluded region changes rather drastically with 
ioff- Although this allows us to conclude that the CoGeNT region is consistent with the XENONIO constraint, the 
region favored by DAMA still appears to be favored by XENONIO. The most important measurements of L^s are 
in the energy bins near XENONlO's threshold, since these determine the lowest recoil energies that XENONIO can 
probe. For illustration, we also show in Fig. [5] the constraint taking the —^Lcs values only in the three bins centered 
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on 4.9, 5.7, and 6.4 keV (and —jLcS at 3.9 kcV so that scintillation always increases with recoil energy), and central 
values elsewhere. 

The second significant effect that we have taken into account is shifts from uncertainty in the halo model. It is 
generally considered that vq = 220km/s is a useful starting point for the halo model; however, if the distribution 
of velocities is broader, then lighter dark matter will be favored, bringing the DAMA region below the XENONIO 
constraint curve. Taking vq = 270 km/s, we see in the right frame of Fig. [2] that the region consistent with DAMA 
shifts to lower masses, and becomes consistent with both the constraints of the null experiments and with the CoGeNT 
preferred region. 
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FIG. 2: Left panel: To illustrate the effect of uncertainties in the XENONIO scintillation factor, LcS, we show the CoGeNT 
and DAMA allowed regions alongside constraints from the new XENONIO analysis using L^a = 0.19 (light gray dashed), and 
I/e// at the central and lower la values from the new measurement [S^ (the two thick black dashed curves). In between the 
central and lower la L^ff curves we have taken — |cr -Le// values in only the three energy bins near XENONlO's threshold 
(black dotted). We also show constraints from the CDMS silicon analysis (red dot-dashed), and the region in which the two 
events observed by CDMS can be produced (blue long dashed). Here, we have used no — 220km/s and Uosc = 500km/s. With 
the lower values of Left values, the tension between XENONIO and CoGeNT is alleviated. Right panel: To illustrate the effect 
from changing the halo model, CoGeNT, CDMS, XENONIO, and DAMA results are shown, but with vo = 270km/s and 
=490km/s. 



To further investigate the source of the remaining tension, we examine the effects of both removing the lowest bin 
from the DAMA data sample, and making the errors in those bins a factor of two larger. In the left frame of Fig. [31 
we have modified the DAMA constraint by showing the region when lowest bin is neglected altogether (grey, 90% and 
99% constraints), as well as the region when the errors are doubled on the lowest two bins (purple, blue, for 90% and 
99%, respectively). We can see that most of the remaining tension is removed: much of the inconsistency seems to 
stem from the lowest one to two bins in the DAMA data. 

In the results shown so far, we have assumed that the couplings of the dark matter to protons are equal to those to 
neutrons. If the dark matter couples to neutrons more strongly than to protons, we can further relax the constraints 
of the CDMS silicon analysis. For example, taking /„ = 1, /p = wc find the results shown in the right frame of 
Fig.H 



Experiment 


Target 


Exposure (kg-d) 


Threshold 


Ref 


CDMSTI 


Ge 
Si 


612. 
83.3 


10 keV 
7 keV 


[33] 
[15, i6] 


XENONIO 


Xe 


316 


2 keV with L^// = 0.19 


[1^ 



TABLE I: Relevant features of the null experiments used in our analysis. Efficiencies and cut acceptances are taken from the 
references above. 
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FIG. 3: Left panel: The DAMA region and XENONIO constraint are the same as in the right panel of Fig.[2]but, for iUustration, 
with the effects of removing the lowest bin of data from the DAMA region (gray), or with the effects of increasing the errors in 
those bins by a factor of two (purple and blue for 90% and 99%). We see that the DAMA region is shifted to smaller masses 
and becomes more consistent with the XENONIO constraint. Right panel: The same in in the left frame, but with dark matter 
couplings only to neutrons, fn = l and fp = 0. The constraints from CDMS silicon relative to the DAMA region are weaker 
than for /„ = fp. Only the 90% confidence limit for DAMA removing the lowest bin is shown. 



III. MODELS OF LOW MASS DARK MATTER 

We begin this section by carrying out a relatively model independent operator analysis of the low mass window 
relevant for CoGeNT, DAMA and CDMS. We then turn to discussing two classes of models: those arising from models 
of Asymmetric Dark Matter (ADM), and from supersymmetry.'^ 

A. Operator Analysis 

The operators connecting dark sector Dirac fermions x or complex scalars <f> to visible sector fermions / are 



Cfs = ^XXff (7) 

Cfv = ^xrxHJ (8) 

C.s = ^'^*<^// (9) 

Csv = ^r^M'/'/V/- (10) 



^ Singlet scalars [34j . Mirror models [35| | and "WIMPless" models [36l | also give rise to light WIMPs, though we do not discuss those 
models here. We refer the reader to these references for details. 
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Fermionic dark matter 



For fermionic dark matter, the scattering cross-section obtained from these operators is given by: 

4a 



'DM'"- AT 



TT (mDM + rriNf 



[Zfp + {A- Z)UY , 



(11) 



where A and Z arc the atomic mass and atomic number of the target nuclei, and a is a number dependent on whether 
the fcrmion is Dirac of Majorana. For Majorana fermions a = 1, while for Dirac fermions, a = 1/4. The effective 
couplings to protons and neutrons, /p,n, can be written in terms of the WIMP's couplings to quarks, and depend both 
on the spin of the DM particle and the mediator. In terms of the dark matter's effective couplings, we can write 



fp,' 



q—u,d.s 



(12) 



q—c.h,t 



where Gq denotes the dark matter's effective coupling to a given quark species. The first term reflects scattering with 
light quarks, while the second term accounts for interactions with gluons through a heavy quark loop. The values 
of f^'"^ are proportional to the matrix element, (qq), of quarks in a nucleon. In our numerical calculations, we use 
values for these quantities based on recent lattice QCD results [13, HI]. If the dark matter scatters with nuclei through 
a scalar interaction, both of the terms in Eq. [T^] typically yield sizable contributions to the cross section. If, on the 
other hand, the dark matter scatters through a vector interaction, then heavy-quark loop contribution is negligible, 
and the nucleon level couplings simplify to 



f _ c)Gu Gd 



f _ Gu Gd 



(13) 



The effective operator relates to the parameters in the underlying Lagrangian through 

Gf _ gDM9f 



(14) 



where ip denotes the mediator, and 5dm, 3/ denote the mediator's couplings to the dark matter and standard model 
fcrmion, respectively. 

The annihilation cross-sections for Dirac fermions XX ~^ f f ■• through scalar, pseudoscalar and vector mediators, 
are 



^ ann 



(15) 



fp 

(7 V 

^ ann 



fV 

^ ann ^ 



An 



f 



(1 - mj/ml,]^j) 



1/2 _ 



''DM 

/ 



1 - m'j/2mlj^j 



G}cf{l~m}/ml,,y/^ 



2 + 



^ 

''DM 



+ 5—1^ 



24(1 



(16) 



(17) 



where c/ is a color factor. 



2. Scalar Dark Matter 



There are related expressions for scalar dark matter scattering and annihilation cross-sections: 

™2 



a 

a — — 



I'M 



TT (™DM + TOTV)' 



(z/p + {A- Z)UY 



(18) 



with 



/p,n — ^ ^^/: 

q—u,d.s 



-^q r{p.n) '^p.n 2 „(^p^n) ^q '^p,n 



(19) 
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for scalar mediators, and 

for vector mediators, a = 1 for real scalars, and a = 1/4 for complex scalars. 

For scalar mediators, the effective coupling is related to the underlying parameters in the Lagrangian through 

Ff _ gumyfiH) 

Ml ' ^'^^ 

where v is the Higgs vev, vj\f2 = {H) = 174 GeV, and Mh is its mass, derived from a term in the Lagrangian 

C = gj^uH^HS^S + VfHff. (22) 

For vector mediators (and complex scalars), the effective coupling is related to the underlying parameters through 

Ff _ goMff/mpM , . 

Ml ■ ^'^^ 

The corresponding annihilation cross-sections for a complex scalar cfxj) — > //, through real scalar and vector medi- 
ators, are 



1) 



(24) 



= ^ E ^/-/(l - ( 1 + ] . (25) 
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3. CoGeNT and Thermal WIMPs 



Beginning with the case of a fermionic WIMP, we can consider either scalar or vector interactions with quarks to 
produce a spin- independent elastic scattering with nuclei. Following Ref. [s^, we show results for each of these cases 
in Fig. 21 In the upper left frame, we have assumed that the dark matter's effective scalar couplings to fermions are 
proportional to the mass of the fermion, Gf oc uif, but are otherwise universal. To produce an elastic scattering 
cross section capable of generating the excesses observed by CoGcNT and DAMA (shown as a dotted region), we 
need an effective coupling of Ges x (1 GcV /mf) ^ 2 x 10~^ GcV~^. If this is the only interaction experienced by the 
WIMP, however, it will be significantly overproduced in the early universe, resulting in a thermal relic abundance in 
excess of the measured abundance of dark matter. The solid lines in each frame of Fig. [5] denote the effective coupling 
strength that leads to the desired thermal dark matter abundance, assuming no other interactions further deplete the 
relic density. For details pertaining to the dark matter's annihilation cross section and the thermal relic abundance 
calculation, see Ref. f39| . 

This apparent conflict between the couplings required to produce the desired elastic scattering rate and the measured 
dark matter abundance can easily be resolved, however, if the fermionic WIMP also annihilates through other process, 
such as through pscudoscalar or axial interactions (which do not, in turn, generate a contribution to the spin- 
independent elastic scattering cross section). In the upper right frame of Fig. 21 we show the effective pseudoscalar 
coupling that leads to the desired dark matter abundance. A 9 GeV WIMP, for example, with effective scalar 
and pseudoscalar couplings of Ges ^ 3 x 10^^ GeV~^(m//l GeV) would satisfy both the elastic scattering and relic 
density requirements. Alternatively, one could also consider effective couplings that are not universal across all fermion 
species. If the dark matter couples preferentially to leptons, for example, its relic abundance could be reduced without 
increasing the clastic scattering cross section. 

As a simple example, we could consider a light fermionic dark matter particle that couples to quarks through the 
exchange of the Standard Model Higgs boson. This scenario, however, is unlikely to be able to generate the excesses 
reported by CoGeNT and DAMA. In particular, even if we take the dark matter's couplings to the Higgs to be of 
order unity, the effectively couplings to quarks would only be GEff ^ 4 x 10"'' GeV~^(m//l GeV), or approximately 
an order of magnitude smaller than is needed to generate the observed number of events. From this example, we see 
that the mediator must either be considerably lighter than the Standard Model Higgs collider limits, or must have 
coupling to quarks that are larger than those of the Standard Model Higgs. 
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FIG. 4: The effective dark matter-fermion couplings required to generate a given spin-independent elastic scattering cross 
section with nucleons (dashed lines), and to produce a thermal abundance consistent with the observed dark matter density 
(solid lines). In each frame, the dotted region represents the approximate range of couplings and masses that leads to a 
signal consistent with that reported by CoGeNT. The various frames show results for a fermionic dark matter particle with 
scalar interactions (upper left), pseudoscalar interactions (upper right), vector interactions (center), or a scalar dark matter 
particle with scalar interactions (lower left), or vector interactions (lower right). In the case of a scalar dark matter particle 
with approximately universal scalar couplings to Standard Model fermions, the coupling size required to generate the excesses 
observed by CoGeNT and DAMA is also approximately the size required to thermally generate the observed density of dark 
matter in the early universe. 



If the dark matter is a Dirae fermion, vector eouplings can also generate a spin-independent elastic scattering cross 
section. In the center frame of Fig. 21 we show the effective couplings required in this case. Here, we have assumed 
that Ges is the same for all species of Standard Model fermions. Once again, we find that the couplings required to 
generate the desired elastic scattering rate are not sufficient to yield the observed dark matter abundance, and other 
interactions are required to avoid the overproduction of dark matter in the early universe. 

If instead of a fermion, the WIMP is a scalar particle, we again can rely on either scalar or vector interactions as 
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the source of the CoGeNT and DAMA excesses. In the case of scalar interactions, we find that the effective coupling 
required to generate the elastic scattering cross section implied by CoGeNT and DAMA also leads to a thermal relic 
abundance consistent with the observed dark matter density (see the lower left frame of Fig. HJ. So, in this simplest 
case we have considered (a scalar dark matter particle with approximately universal scalar couplings to Standard 
Model fermions), we find that if we set the couplings to the size required to generate the observed density of dark 
matter, then we also predict an elastic scattering rate that can easily account for the CoGeNT and DAMA excesses. 
A scalar WIMP with vector interactions, on the other hand, requires another interaction to avoid being overproduced 
in the early universe (see the lower right frame of Fig. U]). 



B. Asymmetric Dark Matter 



In most dark matter models, the abundance of dark matter and the baryon asymmetry of the universe are unrelated 
to each other. The dark matter density is generally determined by its mass and self-annihilation cross section, which 
determine the temperature at which is thermally freezes out in the early universe. The baryon asymmetry, in contrast, 
is set by CP violating phases, and by the strength of the electroweak phase transition. Within this paradigm, the 
densities of dark matter and of baryons in the universe have nothing to do with one another. But despite this 
expectation, observations have revealed that the cosmological densities of these two sectors are similar: 

^ « 5. (26) 
Pb 

Within the context of most thermal dark matter models, one must view this similarity as a coincidence (unless 
anthropic arguments are invoked). 

Motivated by this a ppa rent coincidence, models in which the dark matter abundance is tied to the baryon asymmetry 
have been developed jl2l [l3|. In ADM models, a mechanism (in some cases a higher dimension operator) enforces a 
relationship between the dark matter and baryon chemical potentials: 

ci{nx ~ Tlx) ^ rib - ni, (27) 

where ci is a number 0(1) whose precise value depends on the operator transferring the asymmetry. 

Once the symmetric thermal abundance, nx + nx > has annihilated away, we are left with a relationship between 
the asymmetric relic dark matter density and the relic baryon number density, implying 

moM = SciTOp, (28) 

which combined with Eq. ((26|) requires that the dark matter's mass is roughly within the range of ^1-10 GeV. 

In Ref. [l^l, a class of models was discussed in detail that fits naturally within the CoGeNT, DAMA and CDMS 
windows. In one such model (a supersymmetric example), the dark matter and baryon chemical potentials are related 
to each other through a superpotential interaction 

This operator transfers a lepton asymmetry to the dark matter and predicts a dark matter mass of approximately 
mx ^11 GeV. Other classes of models have also been discussed in which the dark matter is related to the baryon 
density through the out-of-equilibrium decay of a new heavy particle to both sectors. 

ADM thus predicts dark matter with a mass near or within the window implied by CoGeNT and DAMA. In order to 
provide an explanation for these excesses, however, such a dark matter candidate must also possess an appropriately 
large elastic scattering cross section. The magnitude of this cross section is determined not only by the properties 
of the dark matter particle itself, however, but also on the mass and couplings of the mediator, which makes it a 
somewhat more model dependent quantity. As one possibility (as discussed in Ref. [13 )> a simple Z' mediator with 
mass around the TeV scale gives rise to a scattering cross section of the magnitude observed by CoGcNT and DAMA. 
Alternatively, a singlet scalar that mixes with the Higgs may also give rise to an appropriately sized elastic scattering 
cross section. This singlet may be responsible for the mass of the singlet dark matter through an operator 

Wm = XSXX, (30) 

and a mixing with the Higgs through the usual NMSSM term 



(31) 
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gives rise to the appropriate size scattering cross section for CoGcNT and DAMA p^ : 

. 3xlO-cn.^f^Vf-VfAVf^Vf^^Vf^V^ (32) 
where Nn arises from the effective couphng of the Higgs to the target nucleus.^ 



C. (Non-Minimal) Neutralino Dark Matter 



Within the context of the minimal supersymmetric standard model (MSSM), light neutralino dark matter can 
coherently scatter with nuclei through scalar Higgs exchange (and to a lesser extent, through squark exchange), while 
annihilating fairly efficiently through pscudoscalar Higgs exchange (which, unlike scalar exchange, is not s-wave sup- 
pressed). The composition and couplings of a neutralino lighter than mz/2 is constrained by the LEP measurements 
of the invisible Z branching fraction, which require rz^x°x° ^ ^ MeV (at 2a C.L.) (40| . This approximately trans- 
lates to a bound of jiVfs ~ -^mI ^ 0-13 at 2a, where Nf^ and Nf^^ are the up-type and down-type higgsino fractions 
of the lightest neutralino. If we maximize the neutralino-ncutralino-Higgs couplings by saturating this bound, and 
impose the constraints from the Tevatron for Higgs boson production at large tan /3, it has been shown that the largest 
scattering cross-section that can be achieved is ~ 2 x 10""'^ cm^ [i^l. This is near to the CoGeNT and DAMA regions 
[IS SI] J but certainly on the lower edge. 

To obtain larger effective couplings (and thus larger elastic scattering and annihilation cross sections), one could 
consider interactions mediated by Higgs bosons which are lighter than those allowed within the MSSM. Within the 
context of the next-to-minimal supersymmetric standard model (NMSSM), for example, light (~1-100 GeV) scalar 
and pseudoscalar Higgs bosons can exist, while not conflicting with existing collider constraints, as a result of their 
reduced couplings to Standard Model fields. In particular, if the lightest scalar and pseudoscalar Higgs bosons 
contain a significant admixture of the a Higgs singlet, such particles can be much lighter than the MSSM Higgs limits, 
increasing the corresponding effective couplings of the lightest neutralino psl l46j . 

As an example, we could consider a 5-10 GeV neutralino which is mostly bino, but with a ^5% higgsino admixture. 
If such a particle annihilates through the exchange of a < 70 GeV mixed singlet-MSSM-like pseudoscalar Higgs, and 
if tan/3 is relatively large (~30-50), then we find that this process can yield the desired relic abundance. Furthermore, 
if a down- type scalar Higgs (i?i) were also of a comparable mass, and comparably mixed with the Higgs singlet, it 
would mediate an elastic scattering cross section consistent with the rate observed by CoGeNT and DAMA. This 
example serves as an illustration of the more general conclusion that the range of masses and scattering cross sections 
implied by CoGeNT and DAMA requires a particle or particles that are either relatively light or possess relatively 
large couplings to mediate the interactions of the dark matter. 



IV. IMPLICATIONS FOR INDIRECT DETECTION 



Regions with high densities of dark matter are in many cases predicted to produce potentially observable fluxes 
of dark matter annihilation products. As the annihilation rate of dark matter scales with the square of the inverse 
mass of the particle, the prospects for indirect detection are particularly encouraging for dark matter masses as light 
as we are considering in this paper. In this section, we discuss some of the implications of the CoGeNT and DAMA 
excesses for indirect dark matter searches using neutrinos, gamma rays, and charged cosmic rays. We also comment 
on the possibility that white dwarf stars could be used as a probe of dark matter in models capable of generating the 
observed excesses. 



* Note that the reUc abundance arguments described in Sec. IIII Al do not apply in the case of asymmetric dark matter. Furthermore, the 
indirect detection prospects described in Sec. llVl do not apply to this class of dark matter candidates. 
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A. Neutrinos From Dark Matter Annihilations In the Sun 



If dark matter particles are present in the vicinity of the solar system, they will scatter elastically with and become 
captured in the Sun at a rate approximately given by [47| 

- - ^ (ofe) (^) (^) 

where po is the local dark matter density, v is the local root-mean-square velocity of halo dark matter particles, and 
todm is the dark matter mass. The sum is over mass number of the chemical species present in the Sun, and the 
quantity Fa contains information pertaining to the solar abundances, dynamical factors, and form factor suppression 
of each clement {Fn = 1, i^uc ~ l-l; etc.). The quantity 5 is a kinetic suppression factor given by 



where 



Aix) 



3/2 



2/3 



(34) 



The evolution of the number of dark matter particles in the Sun, N, is described by 

TV = C® - A® _ ^0^^ (3g) 

where C® is the capture rate, A® is the annihilation cross section times the relative WIMP velocity per unit volume, 
and i5® is the inverse time for a WIMP to exit the Sun via evaporation (which is negligible for dark matter particles 
heavier than approximately 3-4 GeV 0,14^). A® can be approximated by 

A- = ^, (37) 

where Vcs is the effective volume of the core of the Sun determined roughly by matching the core temperature with 
the gravitational potential energy of a single WIMP at the core radius. This was found in Refs. [H [sq] to be 

%s = 5.7 X 10^0 cm^ ( i^V^' . (38) 

V '"DM / 

Neglecting evaporation, the present WIMP annihilation rate in the Sun is given by 

F ^A^N^ ^l*^^ ^^^^^ (i0VC®A0) , (39) 

where « 4.5 billion years is the age of the solar system. For a 5-10 GeV WIMP, capture-annihilation equilibrium 
is reached so long as the WIMP's annihilation cross section is larger than av > 10^'^*' cm'^/s. Once equilibrium is 
reached, the final annihilation rate (and corresponding neutrino fiux and event rate) has no further dependence on 
the dark matter particle's annihilation cross section. 

As the dark matter particles annihilate, they can generate neutrinos through a wide range of channels. Annihilations 
to bottom quarks, charm quarks, and tau leptons each generate neutrinos in their subsequent decays. In some models, 
dark matter particles can also annihilate directly to neutrino pairs. Once produced, neutrinos travel to the Earth where 
they can be detected. For the mass range of interest here, the strongest constraints on the capture and annihilation 
rate in the Sun come from the Super-Kamiokande experiment|5l| . In Fig. [SJ we show the constraints on a light dark 
matter particle from Super-Kamiokande (adapted from Ref. (S^). From this figure, we conclude that if a dark matter 
particle responsible for the excess observed by CoGeNT annihilates primarily to neutrinos or tau leptons, it should 
have produced a fiux of neutrinos observable to Super Kamiokande. If the dark matter annihilates significantly to 
other types of particles, however, it could produce the rates observed at CoGeNT and DAMA without coming into 
conflict with Super Kamiokande's observations. Alternatively, if the dark matter's annihilation cross section is highly 
suppressed (such as in the case of asymmetric dark matter, as discussed in Sec. IIII B[) . then these constraints do not 
apply. 

Note that the results shown in Fig. [5] differ somewhat from those found in Ref. [S^]. This is because, in Ref. j53 |. 
only scattering off of hydrogen and helium nuclei was included, and a simplified treatment of the velocity distribution 
was adopted [53|. The capture rates used to obtain the limits found in Fig. [5] are in good agreement with the publicly 
avafiable code DarkSUSY lH. 
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FIG. 5: The upper limit on a light WIMP's spin-independent elastic scattering cross section with nuclei from the Super 
Kamiokande experiment, for various choices of the dominant annihilation channel. Also shown is the region of the plane which 
provides a good fit to the excess observed by CoGeNT. From this figure, we conclude that Super Kamiokande excludes a dark 
matter particle which annihilates entirely to neutrinos or taus as a source of the CoGeNT excess. This figure was adapted from 
Ref. [Hi]. 



B. Gamma Ray Searches From Dark Matter 

In generality, the flux of gamma rays from dark matter annihilations can be written as 

*,(E„^,.l<,|„|^^_i^£,^., (40) 

where, a\v\ is the dark matter's annihilation cross section (times relative velocity), ip is the angle observed, p is the 
dark matter density, and the integral is performed over the line-of-sight being observed. dN^/dE^ is the gamma 
ray spectrum generated per WIMP annihilation, which depends on the mass of the dark matter particle and on its 
dominant annihilation channels. The integral in Eq. |40] depends on the distribution of dark matter, and is largest 
in the directions of the sky where very dense and relatively nearby concentrations of dark matter are present. The 
brightest source of gamma rays from dark matter annihilation is generally expected to be the region around the center 
of the Milky Way [s^ . Nearby dwarf spheroidal galaxies are also very promising sources for indirect detection [5^ . 

To estimate the gamma ray flux from the Galactic Center region, we adopt a dark matter distribution which follows 
the commonly used Navarro-Frenk- White (NFW) profile [57l |: 

'^'^={rlR)[lUrlR)f^ (41) 

where r is the distance from the Galactic Center, i? ~ 20 kpc is the scale radius and po is fixed by imposing that the 
dark matter density at the distance of the Sun from the Galactic Center is equal to the value approximately inferred 
by observations (~0.3 GeV/cm"^). 

In Fig.ini we show the spectrum of gamma rays from dark matter annihilations in the inner Milky Way, and compare 
this to the spectrum measured by the Fermi Gamma Ray Space Telescope (FGST) Q. The curves shown in each 
frame are for the case of a 9 GeV dark matter particle annihilating with a cross section of cru = 3 x 10~^^ cm^/s (the 
value estimated for a typical thermal relic) to either hh or t+t~. The left and right frames show the flux predicted 
from within the inner 0.5° and 3° around the Galactic Center, respectively. The error bars shown correspond to the 
measurements of FGST over its first 14 months of observation (as shown in Ref. Q). Each curve has been corrected 
to account for leakage due to FGST's finite point spread function (see Ref. [§1 for details). 

While the fiuxes shown in Fig. |6] do not exceed those measured by FGST, they are quite close ~ within a factor 
of 2-3 below 1-2 GeV. As more data is accumulated, and as the astrophysical backgrounds in this region become 
increasingly understood, FGST will likely become sensitive to 5-10 GeV dark matter particles with annihilation cross 
sections near that predicted for a simple thermal relic. The lone caveat to this conclusion is that, if the dark matter 
profile describing the inner galaxy does not contain a significant cusp (such as that described by the NFW model), 
then the gamma ray fiuxes from dark matter annihilations could be considerably suppressed. 



14 




0.5 1.0 5.0 10.0 50.0100.0 0.5 1.0 5.0 10.0 50.0100.0 

(GeV) (GeV) 

FIG. 6: The spectrum of gamma rays from 9 GeV dark matter particles annihilating in the inner Milky Way, compared to 
the measurements of the Fermi Gamma Ray Space Telescope (FGST). Here, we have assumed a dark matter distribution that 
follows an NFW profile, and an annihilation cross section of av = 3 x lO"'^^ cm'^/s. See text for more details. 

Recently, the FGST collaboration has reported its first limits on the gamma ray flux from a number of dwarf 
spheroidal galaxies, and used this information to constrain the properties of dark matter js^l- Although results are 
not presented for dark matter particles as light as 5-10 GeV (except in the special case of annihilations to 
the FGST collaboration find that dark matter particles with a 30 GeV mass and which annihilate to bb are restricted 
to a\v\ < 6 X 10^^^ cm'^/s, assuming an NFW profile shape. Furthermore, the annihilation products of 5-10 GeV 
dark matter particles are predicted to contribute to the diffuse isotropic gamma ray flux at a level near the current 
sensitivity of FGST 

C. Indirect Detection With Charged Cosmic Rays 

In contrast to most other astrophysical processes, dark matter annihilations are generally predicted to produce 
equal fluxes of matter and antimatter. An excess of antimatter in the cosmic ray spectrum relative to the predictions 
of astrophysical models could thus potentially constitute a signal of dark matter annihilations taking place in the halo 
of the Milky Way. 

Once injected into the interstellar medium, cosmic rays undergo a number of process which can alter their spectra as 
observed at Earth (for a review, see Ref. [fiO])- Such processes include deflection by and diffusion through the Galactic 
Magnetic Field, diffusive reacceleration, convection, and various energy loss processes resulting from interactions with 
gas, radiation fields, and magnetic fields. To confidently interpret the measurements of the cosmic ray spectrum, 
these processes must be understood and adequately modeled. 

At relatively high energies {E > 10 GeV), many of these processes play only a relatively minor role. In particular, 
the effects of convection and reacceleration impact the cosmic ray spectrum much more strongly at low energies. This 
enables one to interpret the high energy measurements of PAMELA and (in the future) AMS-G2 without introducing 
intractable uncertainties associated with the cosmic ray propagation model. For dark matter with a mass in the ~^10 
GeV range, however, the resulting annihilation products will be of sufficiently low energy that they will likely be 
significantly effected by convection and reacceleration, as well as by interactions with the Solar System's magnetic 
field {ie. solar modulation). Such effects may make it difficult to interpret these measurements within the context 
of indirect dark matter searches. That being said, lighter dark matter particles annihilate more often than heavier 
particles, leading to a greater fiux of annihilation products to potentially detect [1^ l6l|. Without exploring these 
competing factors further, we simply comment that light dark matter particles are generally expected to produce 
large fluxes of antimatter in the '^GeV region of the cosmic ray spectrum, but that the complicated nature of cosmic 
ray propagation over this range of energies may make it difficult to clearly identify such a contribution in existing or 
future measurements. 



15 



D. Anomalous White Dwarf Heating 

Dark matter particles with large elastic scattering cross sections (such as those being considered here) can become 
captured by and annihilate in compact stars at very high rates. In relatively cool stars, such as very old white dwarfs, 
the energy injected from these annihilations could potentially compete with or exceed the star's luminosity, leading 
to observable effects [12, HI] . 

As a result of the very high densities of nuclei in white dwarfs, dark matter particles with elastic scattering cross 
sections of cr ~ 10""^*^ cm^ will scatter numerous times as they pass through the volume of the star. In this optically 
thick limit, dark matter particles will be gravitationally captured by a white dwarf at a rate given by 

1 /"^ 

3 / muMV 

where i?wD a-nd A/wd arc the radius and mass of the white dwarf, respectively, pdm is the dark matter density 
in the region of the star, and v is the dark matter's velocity dispersion in that region. For a typical white dwarf 



(i?WD ~ O.OO83i?0, MwD ~ O.7M0). this leads to a contribution to the star's luminosity given by [631 

« 2.5 X lO^^GcV/s 



Pdm \ / 220km/s 



1 GeV/cm^ 



V 



For comparison, the Stefan-Boltzmann law predicts the luminosity of a blackbody to be 

^BB = 47ri?^jjCrT^j-) 

(44) 



Comparing Eqns. l43l and l44l we see that a very old (T ^ 3000 K) white dwarf in a typical region of the galactic disk 
will have only on the order of 0.1% of its luminosity generated by dark matter annihilations. In regions of high dark 
matter density and low dark matter velocity dispersion, however, dark matter annihilations can play a considerably 
greater role. The dark matter density and velocity dispersions within the inner ^10-20 parsecs of dwarf spheroidal 
galaxies is typically on the order of pdm ~ 10^ GeV/cm^ and v ^ 10 km/s, respectively |6J]. In such an environment, 
dark matter of the type required to produce the excesses observed by CoGeNT and DAMA will prevent white dwarf 



stars from cooling below ^ 5000 K |63j . This provides an opportunity for future telescopes not only to confirm a dark 
matter interpretation of the CoGeNT result, but also to potentially map the density and velocity dispersion of dark 
matter in regions such as dwarf spheroidal galaxies. 



V. CONCLUSIONS 



In this paper, we have discussed the recent excess of low energy events reported by the CoGeNT collaboration, 
and considered the possibility that these events, as well as the annual modulation observed by DAMA, are the result 
of a light (ttidm ~ 5 — 10 GeV) dark matter particle. We find that a common dark matter interpretation may be 
compatible with the CoGeNT and DAMA signals, as well as the modest excess recently reported by CDMS. We find 
that this interpretation is also consistent with the null results of XENONIO and the CDMS silicon analysis. 

Dark matter with the properties (mass and elastic scattering cross section) required to produce the CoGeNT and 
DAMA excesses can appear within a variety of theoretical frameworks. From a model independent perspective, 
we note that a ~5-10 GeV scalar dark matter particle with approximately universal scalar couplings to Standard 
Model fermions will naturally produce approximately the observed event rates if its couplings are fixed to obtain the 
desired thermal relic abundance of dark matter. Fermionic dark matter particles are also a viable possibility, although 
multiple types of interactions appear necessary if such a particle is also to satisfy relic density constraints. Within 
the context of specific theoretical models, we have discussed both Asymmetric Dark Matter models and neutralinos 
within extended supersymmetric models, such as the NMSSM. Asymmetric Dark Matter is especially attractive as 
an explanation for the CoGeNT and DAMA signals, as in such scenarios the dark matter's mass is related to the 
proton's mass and is generally predicted to be roughly in the range of muM ~1-10 GeV. 
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We have also discussed the implications of the CoGeNT and DAMA excesses for other dark matter search strategies, 
and find that the prospects for indirect detection are very favorable in many such scenarios. In particular, existing 
limits from low-threshold neutrino detectors, such Super Kamiokande, already constrain the dominant annihilation 
channels of such a WIMP. Searches for dark matter annihilations products with the Fermi Gamma Ray Space Telescope 
(FGST) currently have a level of sensitivity that is very close to being able to test dark matter in a form of a ~5-10 
GeV thermal relic. 



Acknowledgements 

We would like to thank Juan Collar for his assistance in understanding the details of the recent CoGeNT result. We 
also thank B. Feldstcin, A. Pierce, and P. Sorensen for helpful discussions. DH is supported by the US Department 
of Energy, including grant DE-FG02-95ER40896, and by NASA grant NAG5- 10842. ALF and KMZ thank the KITP 
for their hospitality at the workshop "Direct, Indirect and Collider Signals of Dark Matter," where parts of this work 
were initiated. ALF is supported by DOE grant DE-FG02-01ER-40676 and NSF CAREER grant PHY-0645456. 

Bibliography 

[1] R. Bernabei et aZ.. larXiv:1002. 10281 [astro-ph.GA|. 

[2] Z. Ahmed et al. [The CDMS-II Collaboration], 'arXiv:0912. 3592 [astro-ph.CO]. 
[3] C. E. Aalseth et al. [The CoGeNT Collaboration], arXiv: 1002.4703 [astro-ph.CO]^ 

[4] O. Adrian! et al. [PAMELA Collaboration], Nature 458, 607 (2009) arXiviOSlO.lggF [astro-ph]]. 

[5] A. A. Abdo et al. [T he Fermi LAT Collaboration], Phys. Rev. Lett. 102, 181101 (2009) arXiv :0905 .00251 [astro-ph. HE] ] . 
[6] D. P. F inkbeiner, a rXiv:astro-ph/0 409027 ; D. Hooper, D. P. Finkbeiner and G. Dobler, Phy s. Rev. D 76, 083 012 (2007) 
[arXiv:0705. 3655, [astro-ph]]: G. Dobler and D. P. Finkbeiner, Astrop hys. J. 680, 1222 (2008) [arXiv:0712.T038l [astro-ph]]. 
[7] G. Dobler, D. P. Finkbeiner, I. Cholis, T. R. Slatyer and N. Weiner. larXiv:0910.4583l [astro-ph.HE]. 

[8] L. Goodenough and D. Hooper, 'arXiv:0910.2998 [hep-ph]. 

[9] D. Hooper, P. Blasi and P. D. Serpico, arXiv:0810.1527 [astro-ph]; H. Yuksel, M. D. Kistler and T. Stanev. |arXiv:0810. 27841 
[astro-ph]; S. Profumo, arXiv :08T2T4457l [astro-ph]: F. A. Aharonian, A. M. Atoyan, and H. J. Volk, Astron. Astrophys. 
294, L41 (1995); I. Buesching, O. C. de Jager, M. S. Potgieter, and C. Venter, arXiv:0804.0220 [astro-ph]; L. Zhang, and 
K. S. Cheng, Astron. Astrophys. 368, 1063 (2001). M. Kaphnghat, D. J. Phalen and K. M. Zurek, JCAP 0912, 010 (2009) 
pxXrv:0905.0487 [astro-ph.HE]]; J. P. Harding and K. N. Ab azajian, Phys. Rev. D 81, 023505 (2010) |arXiv:0910.4590l 
[astro-ph.CO]]; C. M. Kelso and D. Hooper, arX iv: 1001. 17491 [astro-ph.HE]. 

[10] For DAMA and positron excesses specifically, see K. M. Zurek, Phys. Rev. D 79, 115002 (2009) arXiv:0811.4429 [hep- 
ph]]. More generally on multi-component dark matter: T. Hur, H. S. Lee and S. Nasri, Phys. Rev. D 77, 015008 (2008) 
[arXiv:0710. 26531 [hep-ph]]. M. Fairbairn and J. Zupan, JCAP 0907, 001 (2009) (arXiv:0810.414 7 [hep-ph]]. S. Profumo, 
K. Sigurdson and L. Ubaldi, JCAP 0912, 016 (2009) arXiv:090 7.4374 [hep-ph]]. I. Cholis and N. Weiner, ;arXiv:0911.4954] 

[11] J. Angle et al. [XENONIO CoUaboration], Phys. Rev. D 80, 115005 (2009) arXiv:0910.3698 [astro-ph.CO]]. 

[12] D. B. Kaplan, Phys. Rev. Lett. 68, 741 (1992); S. Nussinov, Phys. Lett. B 165, 55 (1985); S. M. Barr, Phys. Rev. D 44, 
3062 (1991); S. M. Barr, R. S. Chivukula and E. Farhi, Phys. Lett. B 241, 387 (1990); S. B. Gudnason, C. Kouvaris and 

F. Sannino, Phys. Rev. D 73, 115003 (2006) arXiv:hep-ph/0603014 ; S. Dodelson, B. R. Greene and L. M. Widrow, Nucl. 
Phys. B 372, 467 (1992); M. Fuju and T. Yana gida, Phys. Lett. B 54 2, 80 (2002) arXiv:hep -ph/0206 066 ; R. Kitano 
and I. Low, Phys. Rev. D 71, 02 3510 (2005) arXiv:hep-ph/041 1 133| ; R. Kitano, H. Murayama and M. Ratz, Phys. 
Lett. B 669, 145 (2008 ) 'arXiv:0807.4313 [hep-ph]]; G. R. Farrar and G. Zaharij as, Phys. Rev. Le tt. 96, 041302 (2006) 
arXiv:hep-ph/0510079| . Z. Berezhiani and A. Lepidi, Phys. Lett. B 681, 276 (2009) larXiv:0810."T3T7l [hep- ph]]. G. D. Kribs, 
T. S. Roy, J. Terning and K. M. Zurek, arXiv:0909.2034 [hep-ph]. Y. Cai, M. A. Luty and D. E. Kaplan. larXiv:0909.549"9l 
[hep-ph]. H. An, S. L. Chen, R. N. Mohapatra and Y. Zhang, arXiv:0911.4463 [ hep-ph]. 

[13] D. E. Kaplan, M. A. Luty and K. M. Zurek, Phys. Rev. D 79, 115016 (2009) [arXiv:0901.4TT7l [hep-ph]]. T. Cohen and 
K. M. Zurek, iirXiv:0909. 2035 [hep-ph]. 

[14] J. Angle et al. [XENON CoUaboration], Phys. Rev. Lett. 100, 021303 (2008) larXiv:0706.0039l [astro-ph ]] . 

[15] D. S. Akerib et al. [CDMS Collaboration], Phys. Rev. Lett. 96, 011302 (2006) 'arXiv:astro-ph/0509259" . 

[16] J. Filippino [CDMS Collaboration], Les Rencontres de Physique de la Vallee D'Aost, Nuovo Cimento C 32 05-06 (2009), 
p. 45. 

[17] F. Petriello and K. M. Zurek, JHEP 0809, 047 (2008) [arXiv:08 06.3989^ [hep-ph]]. 

[18] S. Chang, A. Pierce and N. Weiner, Phys. Rev. D 79, 115011 ( 2009) arXiv:0808.0196' [hep-ph]]. 

[19] M. Fairbairn and T. Schwetz, JCAP 0901, 037 (2009) 'arXiv:0808.0704 [hep-ph]]^ 

[20] C. Savage, K. Freese, P. Gondolo and D. Spolyar, JCAP 0909, 036 (200 9) [MXiv:0901. 27131 [astro-ph]]. C. Savage, 

G. Gelmini, P. Gondolo and K. Freese, JCAP 0904, 010 (2009) arXiv:0808.3607" [astro-ph]]. 
[21] E. M. Drobyshevski, Mod. Phys. Lett. A 23, 307;^ ( 2008 ) arXiv:0706.3095 [physics.ins-det]]. 
[22] C. S. Kochanek, Astrophys. J. 457, 228 (1996) [arXiv:astro-ph /9505068 . 

[23] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996) |arXiv:hep-ph/9506380 . 



BIBLIOGRAPHY 



17 



[24] G. Duda, A. Kemper and P. Gondolo, JCAP 0704, 012 (2007) [arXiv:he p-ph/0608035' . 

[25] G. Fricke, C. Bernhardt, K. Heilig, L. A. Schaller, L. Schellenberg, E. B. Shera and C. W. de Jager, Atom. Data Nucl. 
Data Tabl. 60, 177 (1995). 

[26] R. Bernabei et al, Eur. Phys. J. C 53, 205 (2008) [arXiv:0710.0288 : [astro-ph]]. 

[27] A. Manzur, A. Curioni, L. Kastens, D. N. McKinsey, K. Ni and T. Wongiirad. larXiv:0909.1063l [phvsics.ins-det]. 



[28] S. Yellin , Phys. Rev. D 66, 032005 (2002) [arXiv: physics /0203002 
[29] R. Foot, |arXiv :0804.4518 [hep-ph] 
[30] J. Graichen, K. Maier, J. Schuth, A. Siepe and W. von Witsch, Nucl. lustrum. Meth. A 485, 774 (2002) [Erratum-ibid. A 

501, 630 (2003)]. 

[31] J. Kopp, T. Schwetz and J. Zupan, JCAP 1002, 014 (2010) [arXiv:0912.4264l [hep-ph]] . 

[32] P. Sorenson, PhD Thesis. 

[33] Z. Ahmed et al. [CDMS Collaboration]. l arXiv:0802. 35301 [astro-ph]. 

[34] C. P. Burgess, M. Pospelov and T. ter Veldhuis, Nucl. Phys. B 619, 709 (2001) 'arXiv: hep-ph/00 11335'. S. Andreas, 
T. Hambye and M. H. G. Tyt gat, JCAP 0810, 034 (2008) arXiv:0808.0255 [hep-ph]]. Y. G. Kim and S. Shin, JHEP 0905, 
036 (2009) |arXiv:0901.2"609l [hep-ph]] . ' 

[35] R. Foot, Int. J. Mod. Phys. D 13, 2161 (2004) |arXiv:astro-ph/0407623] . R. Foot, Phys. Rev. D 78, 043529 (2008) 



|arXiv:0804.4"5T8l [hep-ph]] 

[36] J. L. Feng and J. Kumar, Phys. Rev. Lett. 101, 231301 (2008) arXiv :0803.4T96l [hep-ph]]. J. L. Feng, J. Kumar and 

L. E. Strigari, Phys. Lett. B 670, 37 (2008) arXiv:0806.3746 [hep-ph]]. 
[37] J. Giedt, A. W. Thomas and R. D. Young, arXiv:0907.4177 [hep-ph]^ 

[38] R. D. Young and A. W. Thomas, Phys. Rev. D 81, 014503 (2 010) ^X iv:0901.3310'[hep-lat]]; D. Toussaint and W. Freeman 
[MILC Collaboration], Phys. Rev. Lett. 103, 122002 (2009) [arXiv:0 905.2432 [hep-lat]]. 



[39] M. Behran, D. Hooper, E. W. Kolb and Z. C. Krusberg, Phys. Rev. D 80, 043509 (2009) arXiv:0808.3384 ' [hep-ph]]. 
[40] C. Amsler et al. (Particle Data Group), Phys. Lett. B 667, 1 (2008). 

[41] V. M. Abazov et al. [DO Collaboration], arXiv:0912.0968 [hep-ex]; V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 
102, 051804 (200 9) arXiv:08 11.0024 [hep-ex]]; V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 10 1, 221802 (2008) 
|arXiv:0805.3556l [hep-ex]]: V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 101, 071804 (2008) |arXiv:0805.249n 
[hep-ex]]; T. Aahonen et al [CDF Colla boration], arXiv:0906. 1014. [hep-ex]: D. E. Acosta et al. [CDF Collaboration], Phys. 
Rev. D 72, 072004 (2005) [arXiv:hep-ex/0506042|. 

[42] E. Kufiik, A. Pierce and K. M. Zurek, arXiv:1003.068T [hep-ph]. D. Feldman, Z. Liu and P. Nath, arXiv: 1003.0437 [hep-ph]. 

[43] see e.g. A. Bottino, N. Fornengo and S. Scopel, Phys. Rev. D 67, 063519 (2003) arXiv:hep-ph/0212379 . A. Bottino, 
F. Donato, N. Fornengo and S. Scopel, Phys. Rev. D 78, 083520 (2008) arXiv:0806.4099 [hep-ph]]. V. Niro, A. Bottino, 
N. Fornengo and S. Scopel, Phys. Rev. D 80, 095 019 (2009) arXiv:0909.2 348 [hep-ph]]. A. Bottino, F. Donato, N. Fornengo 
and S. Scopel, Phys. Rev. D 72, 083521 (2005) |arXiv:hep-ph/0508270| . 

[44] D. Ho oper and T. Pleh n, "Supersymmetric dark matter: How light can the LSP be?," Phys. Lett. B 562, 18 (2003) 
|arXiv:hep-ph/0212226]; 



[45] J. F. Gunion, D. Hooper and B. McEhath, Phys. Rev. D 73, 015011 (2006) arXiv:hep-ph/0509024 
[46] F. Ferrer, L. M. Krauss and S. Profumo, Phys. Rev. D 74, 115007 (2006) arXiv:hep-ph/0609257 
[47] A. Gould, Astrophys. J. 388, 338 (1991). 
[48] A. Gould, Astrophys. J. 321, 560 (1987). 

[49] K. Griest and D. Seckel, Nucl. Phys. B 283, 681 (1987) [Erratum-ibid. B 296, 1034 (1988)]. 
[50] A. Gould, Astrophys. J. 321, 571 (1987). 

[51] S. Desai et al. [Super-Kamiokande Collaboration], Phys. Rev. D 70, 083523 (2004) [Erratum-ibid. D 70, 109901 (2004)] 

[arXiv:hep-ex/0404025 . 

[52] D. Hooper, F. Petriello, K. M. Zurek and M. Kamionkowski, Phys. Rev. D 79, 015010 (2009) |arXiv:0808."2464l [hep-ph]]. 

[53] J. Edsjo, private communication. 

[54] P. Gondolo, J. Edsjo, P. UUio, L. Bergstrom, M. Schelke and E. A. Bahz, JCA P 0407, 008 (2004) [arXiv :astro-ph/0406204l. 



[55] L. Bergstrom, P. UUio and J. H. Buckley, Astropart. Phys. 9, 137 (1998) arX iv:astro-ph/9712318| ; V. Berezinsky, A. Bot- 
tino and G. Mignola, Phys. Lett. B 325, 136 (1994) arXiv:hep-ph/9402215 



[56] N. W. Evans, F. Ferrer and S. Sarkar, Phys. Rev. D 69, 123501 (2004) arXiv: ast ro-ph/03 1 1 145|| ; L. Bergstrom and 
D. Hooper, Phys. Rev. D 73, 063510 (2006) arXiv:hep-ph/0512317l; G. D. Martinez, J. S. Bullock, M. Kaplinghat, 
L. E. Strigari and R. Trotta, JCAP 0906, 014 (2009) arXiv:0902.47T5' [astro-ph. HE]]. 

[57] J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. J. 462, 563 (1996) [arXiv:astro-ph /9508025] ; J. F. Navarro, 
C. S. Frenk and S. D. M. Wh ite, Astrophys. J. 490, 493 (1997). 

[58] A. A. Abdo et a/., larXiv: 1001 .4531 [astro-ph. CO]. 

[59] K. N. Abazajian, P. Agrawal, Z. Chacko and C. Kihc, arXiv: 1002:3820" [astro-ph.HE]. 



[60] A. W. Strong, I. V. Moskalenko and V. S. Ptuskin, Ann. Rev. Nucl. Part. Sci. 57 (2007) 285 la rXiv:astro-ph/0701517 . 
[61] A. Bottino, F. Donato, N. Fornengo and S. Scopel, Phys. Rev. D 77, 015002 (2008) arXiv:0710.0553l[hep-phJj: A. Bottino, 

F. Donato, N. Fornengo and P. Salati, Phys. Rev. D 72, 083518 (2005) |arXiv:hep-ph/0507086i ; 
[62] I. V. Moskalenko and L. L. Wai, Astrophys. J. 659, L29 (2007) [ar Xiv:astro-ph/0702654 |; G. Bertone and M. F airbairn, 

Phys. Rev. D 77, 043515 (2008) arXiv:0709.1485 [astro-ph]]; M. McCuUough and M. Fairbairn, larXTv:1001. 27371 [hep-ph]. 

[63] D. Hooper, D. Spolyar, A. Vallinotto and N. Y. Gnedin, arXiv:1002.0005 [hep-ph]. 

[64] M. G. Walker, M. Mateo, E. W. Olszewski, O. Y. Gnedin, X. Wang, B. Sen and M. Woodroofe. rarXiv:0708.00T0l [astro-ph] . 



